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Abstract. The eigenvalue density for members of the Gaussian orthogonal and unitary en- 
sembles follows the Wigner semi-circle law. If the Gaussian entries are all shifted by a constant 
amount c/(2Af) 1 / 2 , where N is the size of the matrix, in the large N limit a single eigenvalue 
will separate from the support of the Wigner semi-circle provided c > 1. In this study, using 
an asymptotic analysis of the secular equation for the eigenvalue condition, we compare this 
effect to analogous effects occurring in general variance Wishart matrices and matrices from 
the shifted mean chiral ensemble. We undertake an analogous comparative study of eigenvalue 
separation properties when the size of the matrices are fixed and c — » oo, and higher rank ana- 
logues of this setting. This is done using exact expressions for eigenvalue probability densities 
in terms of generalized hypergeometric functions, and using the interpretation of the latter 
as a Green function in the Dyson Brownian motion model. For the shifted mean Gaussian 
unitary ensemble and its analogues an alternative approach is to use exact expressions for 
the correlation functions in terms of classical orthogonal polynomials and associated multiple 
generalizations. By using these exact expressions to compute and plot the eigenvalue density, 
illustrations of the various eigenvalue separation effects arc obtained. 



1. Introduction 

1.1. The aims of the paper. Our interest is in the shifted mean Gaussian ensemble, the general 
variance Laguerre ensemble, and the shifted mean chiral ensembles from random matrix theory. 
The shifted mean Gaussian ensemble consists of real symmetric (/3 = 1) or complex Hermitian 
{(3 = 2) matrices H with joint distribution of the elements proportional to 

(1-1) exp(-^Tr(ff-i^)) 2 ) 

where H®> is a fixed matrix of the same type as H . It follows immediately from (|1 . ip that 
H = G + where G is a member of the (zero mean) Gaussian orthogonal ensemble (/? — 1) 
or Gaussian unitary ensemble {(3 — 2). In the early stages of the development of random matrix 
theory in nuclear physics, Porter studied the real case of Ijl.ip . with all entries of H^ ' equal 
to a constant /j,, through computer experiments of the eigenvalue distribution. It was found 
that typically all but one of the eigenvalues followed the prediction of the Wigner semi-circle 
law (see e.g. [17]), and so to leading order were supported on [— y2N, \/2N]. The remaining 
eigenvalue — the largest — separated from the semi-circle. We owe our knowledge of this to the 
work of Lang [27] , who followed up on this empirical finding with a theoretical study predicting 
that the largest eigenvalue would occur near fxN . 

Jones et al. [21], by way of their analysis of a spherical model of spin glasses [25], made a 
quantitative study of this effect in the large N limit. They showed (see Section [2] below) that 
with i2"(°) having all elements equal to c/y/2N, the separation of the largest eigenvalue only 
occurs for c > 1, thus exhibiting a phase transition in c. 

Many year later [5], an analogous effect has been exhibited in the general variance Laguerre 
(also known as Wishart) ensemble. This ensemble consists of matrices X^X (Wishart matrices) 
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where the joint distribution on the elements of the n x m (n > m) matrix X is proportional to 
(1.2) exp f - ^Tr^tjQT 1 ) 



(1.5) 



2 

where E is the mxm covariance matrix, and as such is required to have all eigenvalues positive. 
The elements of X are real for (3 = 1, and complex for (3 = 2. Writing Y = XYr 1 ! 2 we see Y^Y 
is a Wishart matrix with E = l m , and furthermore 

(1.3) x*x = y} ,2 y^yy}' 2 ~ r f y£ 

where the tilde denotes the matrices are similar and so have the same eigenvalues. With 

(1.4) E = diag(6, (l)™- 1 ) 

(here the notation (l) m_1 denotes 1 repeated m— 1 times) it was shown that the largest eigenvalue 
separates from the leading support on (0,4m) provided bm > 2. 

It is the objective of this paper to exhibit the mathematical similarities between these two 
eigenvalue separation phenomena. Moreover, we add to this a third example, namely the shifted 
mean chiral matrices 

Omxm X + X^ 

xt + (x(°))t o„ x „ 

where X is an n x m Gaussian random matrix with distribution given by (|1.2p in the case E = l m . 
The matrix X^ is a constant matrix of the same type as X. With all elements of X^ ' equal to 
c/(2 v / m), a separation of the largest eigenvalue in modulus (the eigenvalues of (|1.5|) occur in ± 
pairs) only occurs for c > 1. 

1.2. Related works. In the mathematics literature the problem of the statistical properties of 
the largest eigenvalue of (|1.1|) in the case that all entries of H (°> are equal to [i was first studied 
by Fiiredi and Komlos 20J . In the more general setting of real Wigner matrices (independent 
entries i.i.d. with mean \x and variance a 2 ) the distribution of the largest eigenvalue was identified 
as a Gaussian, so generalizing the result of Lang. Only in recent years did the associated phase 
transition problem, already known to Jones et al. |25j . receive attention in the mathematical 
literature. In the case of the GUE, this was due to Peche [35], while a rigorous study of the GOE 
case can be found in the work of Maida [S3] . 

The paper [5] by Ben Arous, Baik and Peche proved the phase transition property relating to 
(|1.2[) in the complex case with E given by E = diag((6) r , (l) m ~ r ) f (|1.4[) corresponds to r = 1). 
Subsequent studies by Baik and Silverstein .6] , Paul [31] and Bai and Yao [I] considered the real 
case. Significant for the present study is the result of [4], giving that for (|1.2|) with E given as 
above, the separated eigenvalues have the law of the rxr GUE. The case (3 = 4 — not considered 
here — corresponding to self dual quaternion real matrices, is studied in the recent work of Wang 

m- 

The eigenvalue probability density function (|1.1|) is closely related to the Dyson Brownian 
motion model |15] in random matrix theory (see Section 3 below). It is also referred to as a 
Gaussian ensemble with a source. In this context the case of H^ ' having a finite rank has been 
studied by a number of authors [3 [22 El H] ■ However our use of this differs in that we will keep 
N fixed, and exhibit phase separation as a function of the perturbing parameter. 

1.3. Plan. We begin in Section 2 by showing how in each case the criteria for separation can 
be deduced from the eigenvalue equation implied by regarding the parameters b, c, as the pro- 
portionalities in rank 1 (or in the case of (|1.5p . rank 2) perturbations. In Section 3 we discuss 
analytic features of the eigenvalue separation effect in the case that the matrices H^°\ E and X 
are of rank r. Working with matrices of fixed size, we make use of a generalized hypergeometric 
function form of the joint eigenvalue probability density function. 
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In Section 4, in the case (3 = 2 (complex) case, we draw attention to known explicit forms for 
the fc-point correlation functions, and show how these are well suited to exhibiting eigenvalue 
separation both analytically and numerically. The content of Section 5 is to exhibit an analogue 
of eigenvalue separation for the general variance Laguerre ensemble, variance matrices of the 
form (|1.4p . and its generalization £ = diag((fe) r , (l) m ~ r ) in the limit 6^0. An applied setting 
is identified in Section 6 as a motivation for future work. 



2. Secular equations 

Consider first the shifted mean Gaussian ensemble (jl.ip . and for definiteness suppose (3 = 1 
(N x N real symmetric matrices). The simplest case is when all elements of Hq are constant, 
equal to /j, say. Then 

(2.1) H = G + nxx T 

where x is a column vector with all entries equal to 1, and G is a member of the (zero mean) 
Gaussian orthogonal ensemble. Diagonalizing G, G = OLO T , L = diag (a±, . . . , a^v), and using 
the fact that O t x = : y\s distributed as a standard Gaussian vector shows that from the viewpoint 
of the eigenvalues, the right hand side of (|2.ip can be replaced by L + fiyir . We seek the 
eigenvalues of this matrix J3J. 

Lemma 2.1. The eigenvalues of the matrix 

H := diag(ai, . . . , a N ) + /j,yif r 

are given by the zeros of the solution of the secular equation (for the use of this terminology, see 



N 



(2.2) = l-/i^^ 



2 



r-f A - di 

Assuming the ordering a% > • • • > ajy, and that fi > 0, a corollary of this is that the eigenvalues 
satisfy the interlacing 

(2.3) Ai > ai > A 2 > a 2 > • • • > \n > a^v- 
Proof. With A = diag (a i, . . . , ajv) we have 

(2.4) det(ljvA -H) = dot^A - A) det(ljv - nyy T {l N X - ^T 1 )- 
The matrix product in the second determinant has rank 1 and so 



(2.5) det(ljv - ^yy T (l N X - AY 1 ) = 1 - Tr( M yy T (ljvA - Ay 1 ) = 1 - /i 



N 
i=l 



2 



The characteristic polynomial (|2.4[) vanishes at the zeros of this determinant, but not at the 
zeros of det(l n A — A) due to the cancellation with the poles in (|2.5[) . Thus the condition for the 
eigenvalues is the secular equation (|2.2[) . The interlacing condition can be seen by sketching a 
graph, taking into consideration the requirement [iyl > 0. □ 
We are interested in the position of the largest eigenvalue of H or equivalently H as a function 
of fi and N, which unlike the other eigenvalues is not trapped by the eigenvalues of A. Following 
[25], to do this requires using the Wigner semi-circle law for the eigenvalue density of Gaussian 
orthogonal ensemble matrices, 

w f (2N/nJ)(l~xyj 2 ), \x\<J 

(2 - b > 0, \x\>J 
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where J = (27V) 1 / 2 . Thus in (|2.2p we begin by averaging over the components of the eigenvectors 
and thus replacing each yf by its mean value unity. Because the density of {a^} is given by p^, (x), 
for large N (|2.2|) then assumes the form of an integral equation, 

(2.7) = 1- a* / ^^dx. 



-j 



A — x 



We seek a solution of this equation for A > J. 
Now, for A > J, 



,2.x, / ^f-^^^r^w^^ 

- J k=0 J ~ J 



But it is well known (see [H]), and can readily be checked directly, that 
(2.9) J' ' x*p™(x)dx = ^(J/2) 2fc+1 c fc 



where Ck denotes the fc-th Catalan number. Substituting (|2.9|) in (|2.8[) . and using the generating 
function for the Catalan numbers 



oo 1 

^c^ = l(l- (1-4^), 



fc=0 

we deduce the integral evaluation 

(2.10) £*&£*_™( 1 _< 1 _ JW *), W>J. 

Substituting this in (|2.7j) shows a solution with A > J is only possible for /i > J/(2N), and 
furthermore allows the location of the separated eigenvalue to be specified. 

Proposition 2.2. [25] Consider shifted mean Gaussian orthogonal ensemble matrices as specified 
by i lg.il) . and write fj, = c/J, J = (27V) 1 / 2 . For large TV, and with c > 1, a single eigenvalue 
separates from the eigenvalue support in the case \x = 0, (— J, J), and is located at 

(2.11) A= J/ ' 



2 

The case of the shifted mean Gaussian ensemble (jl.ip with (3 — 2 (N x N complex Hermitian 
matrices) and all elements of Hq constant, is very similar. With G a (zero mean) member of the 
Gaussian unitary ensemble, the shifted mean matrices have the same eigenvalues distribution as 
the random rank 1 perturbation 

(2.12) di&gG + fixx^ 

where x is an ri-component standard complex Gaussian vector. Thus with yf replaced by \yi\ 2 
the secular equation (12. 2|) again specifies the eigenvalues. Since the mean of \yi\ 2 is unity, and 
the eigenvalue density of Gaussian unity ensemble matrices is again given by the Wigner semi- 
circle law (|2.6p , all details of the working leading to Proposition 12.21 remain valid, and so the 
proposition itself remains valid for shifted mean unitary ensemble matrices. 

We turn our attention next to correlated Wishart matrices (| 1 ,3|) , and consider for definiteness 
the complex case /? = 2. To begin we make use of the general fact that the non-zero eigenvalues 
of the matrix products AB and BA are equal. This tells us that the non-zero eigenvalues of 
FtyS are the same as those for YTY^ . With E given by (TQ|) . 

(2.13) yEr f = y r f + 
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where y is an n-component standard complex Gaussian vector and Y is the (n — 1) x to matrix 
constructed from Y by deleting the first column. For n > m the eigenvalues of this matrix have 
the same distribution as 

(2.14) diag((0)' 1 - 1 -" 1 , oi, ... , a m ) + byy* 

where {(0) I1_1_m , oi, . . . , a m } are the eigenvalues of YY*, and thus {ai, . . . , a m } are the eigen- 
values of Y^Y. 

We recognize ()2. 14[) as structurally identical to (|2.12p . Hence the condition for the eigenvalues 
is given by (|2.2|) with yf replaced by \yf\, fi replaced by b, N replaced by n — 1, a; = (i — 
1, . . . , n — 1 — m), ai i— > a,_( n _i_ m ) (i = n — m, . . . , n — 1). In the limit that to — ► oo with n—m 
fixed, the density of eigenvalues for complex Wishart matrices follows the Marcenko-Pastur law 
(see e.g. [IT]) 

(2.15) hm ffitym) = ( ^VSDd - < y < 4 
v ; m—oo K W vw ; \ 0, otherwise 

Hence, analogous to ()2.10p , after averaging and for large m with n — m fixed the secular equation 
becomes the integral equation 

(2.16) l = bm -^-dy. 

J a x ~y 

Substituting (|2.15|) . and changing variables y t— > mx 2 reduces this to read 

(2.17) U*/'"-''^'^ 

7r J Q X/m — x z 

Noting that (m/n)(l — .t 2 /4) 1//2 = pWJx) with J = 2, JV = min the latter, following the strategy 
which lead to (|2.10|) the integral can be evaluated to give 

(2.18) l = ^-(l — (1 — 4m/A) 1//2 j. 

This equation only permits a solution for bm > 2, which therefore is the condition for eigenvalue 
separation. Moreover, the position of the separated eigenvalue can be read off according to the 
following result. 

Proposition 2.3. |5j Consider the rank 1 perturbation 12.13]) and suppose b > b/m. In the 

limit vn — ► oo with n — m fixed, provided b > 2 a single eigenvalue separates from the eigenvalue 
support (0, Am) in the case 6 = 0, and occurs at 

b 2 



A = 



6-1 



The same result holds for real Wishart matrices. In this case, in (|2.13[) we must replace yyi by 
yy 1 ^ where y is a standard real Gaussian vector, and we must replace YY> by YY T where YY T 
is a real Gaussian matrix. Only the mean value of yf, and the density of eigenvalues of YY T , 
are relevant to the subsequent analysis and as these are unchanged relative to the complex case 
the conclusion is the same. 

Also of interest in relation to (|2.13p is the limit to — > oo with n/m = 7 > 1. This is distinct 
to the case n — m fixed, because the eigenvalue density is no longer given by (|2.15p , but rather 
the more general Marcenko-Pastur law (see e.g. [17] ) 



(2.19) lim pf^''{x-jm) = ( 1 - -)s(x) H — ^ — —Xc<x<o 
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where 

c=(l-^) 2 , d=(l + V7) 2 . 
Note that the delta function term in (12.19P is in keeping with the fraction of zero eigenvalues of 
Y f f equalling (1 - l/ 7 ) (recall (OP ). 

In this case the averaged large m limit of the secular equation becomes the integral equation 

(2.20) 

A 7 i c X/n-x 

where 

-MP, M ._ y/(x-c)(d-x) 
Proceeding as in the analysis of the integral in (|2.8p . for z :— X/n > d we write 

(2.21) f ^L® = I f) 1 /' ^ <*c. 

7c ^ - x z^z k J c W 

Use can now be made of the known fact (see e.g. [40] ) 

(2-22) jf^'W^^-ElOG*^. 

Z — 1 

where the coefficient of 7* is the Narayana number N(k, i — 1). Dehning 

A fc (p,«) :=^iV(fc,z-l)pV +1 ^ 
i=l 

and introducing the generating function 

00 

i/>(p,q,t) = J2 Ak (p,q)t k 

fe=l 

it is fundamental to the theory of Narayana numbers |38j that 

(2.23) ttp(p, q,t) = ^(l-u-v-(l-2(u + v) + (u- v) 2 ) 1/2 ) 

where u := pt, v := qt. Thus with (|2.22ll substituted in (|2.21jl the sum can be expressed in the 
form (|2.23[) to give 

(2.24) [ ^ dx = i f 1 - 1Z1 _ f ! _ HCT±1) + (t-I^x i/> 
J c z — x 2 V z V 2 z 2 J 

Note that in the case 7 = 1 we reclaim the evaluation of the integral implied by (|2.18l) . 

For z := X/n > d both terms in (|2.20p are decreasing functions of z and so take on their 
maximum value when z = d. The square root in (|2.24[) then vanishes, and we find that for (|2.20[) 
to have a solution we must have 

bn> V7(l + V7)- 

In this circumstance the value of A solving (|2.20[) . and thus the location of the separated eigenvalue 
can be made explicit. 

Proposition 2.4. [5 Consider the rank 1 perturbation 12.13\) and suppose b 1— » b/m. In the 
limit Til — ► 00 with n/m = 7 > 1 fixed, provided b > 1 + 1/^/7 a single eigenvalue separates from 
the eigenvalue support in the case b — 0, and occurs at 

X = b(l ■ 7 " 



6-1 



EIGENVALUE SEPARATION IN SOME RANDOM MATRIX MODELS 



7 



The reasoning that tells us Proposition 12.31 holds in the real case implies too that Proposition 
similarly holds in the real case. 
It remains to study the eigenvalues of (|1.5|) in the case that all entries of X 1 - ' are equal to /i. 
Thus 



(2.25) 



where l p is the p x 1 vector with all entries equal to unity. In general a matrix of the form 
(|1.5[) has n — m zero eigenvalues and 2m eigenvalues given by ± the positive square roots of 
the eigenvalues of Y^Y, Y = X + X^°\ Correspondingly, the matrix of eigenvectors has the 
structured form 



(2.26) 



U n xn 
Vm. x n. 



vm. x vn 



with the eigenvalues ordered from largest to smallest. Here U nxm denotes the matrix U nxn with 
the final n — m columns deleted, while V mxrl denotes V mxm with n — m columns of zeros added. 
The matrix (|2.26[) is real orthogonal when (1 1 . 5[) has real entries, and is unitary when the entries 
of (|1.5j) are complex. 

Substituting (|2~2"5]) for in (TO]) and conjugating by (|2.26[) shows (|1.5p has the same 

eigenvalue distribution as 



(2.27) 



diag(A!,...,A m ,(0r 



-A, 



-A, 



u* 



V n 
-ft 



where here * denotes complex conjugation, and 

tTtt jr 



1 1 TJ 



- \ T v 

— 1 m, v mxm- 



(2.28) 

With the diagonal matrix in (|2.27|) denoted A the characteristic polynomial for the matrix sum 
(|2.27p can be written 



(2.29) \ n - m Y[{\ 2 - Af)det(l 



—>-M 

II,, 



-vt 



(Al n+ro - A)- 1 ). 



Use of the general formula 

det(l p - A pxq B qxp ) = det(l q 
shows that the determinant in (|2.29p is equal to the 2x2 determinant 



Bq X pAp X q) 



(2.30) det (l 2 + fx 

where, with u T = (uj)j=i „ and = ((vj)j= 



an ai2 

0,21 0,22 



an = 



j )j=l,...,mi 
0,22 = 



(o)"-™), 



H x2 - x ) 

3 = 1 3 
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Substituting (|2.30|) in (|2.29[) shows ()1 .5|) again has n — m zero eigenvalues, with the remaining 
eigenvalues, which come in ± pairs, given by the zeros of (|2.30p . From the definitions (|2.27[) and 
§TM) we deduce that 

(«,«,>=<), 2(\v j \ 2 )=2(\u j \ 2 ) = l 

and so after averaging over the eigenvectors, and to leading order in m with n — m fixed, the 
eigenvalue condition reduces to 

m * 

(2.31) l = n^ 



^ A 2 - A 2 ■ 

i=i 3 

The positive solutions of (|2.31[) . {xj}j=i,...,m sav ; exhibit the interlacing 

< Ai < x\ < A 2 < x 2 < ■ ■ ■ < A m < x m 
and we know that the density of the positive eigenvalues is given by (x) with 

(2.32) N i-> to, J = 2y^. 

Hence, after averaging over (|2.31[) becomes the integral equation 

2.33 1 = — / -dy. 

^ Jo (A 2 /4to) - y 2 

This integral is essentially the same as that appearing in (|2.17|) . so the evaluation implied by 
(|2.18p allows the integral in (|2.33p for A 2 /4m > 1 to be similarly evaluated, giving 



1 = ^(i-(1-4to/A 2 ) 1/2 ). 



This is precisely the equation obtained by substituting (|2. L0|) in (|2.7j) . after the identifications 
(|2.32p . Hence we obtain for the shifted mean chiral ensemble a result identical to that obtained 
for the shifted mean Gaussian ensemble given in Proposition [ 



Proposition 2.5. Consider shifted mean chiral ensemble matrices specified by \1.5\) , and set 
/i = c/J, J := 2y / ro. For large m, with n — m fixed and c > 1, a single eigenvalue separates from 
the support (0, J) of the positive eigenvalues, and its location is specified by k2.11\) . 



3. Green function viewpoint 

The shifted mean Gaussian and shifted mean chiral random matrices can be used to formulate 
the Dyson Brownian motion of random matrices (see e.g. [17]). As we will see, this viewpoint 
can be used to exhibit the separation of a set of r eigenvalues, in the setting that the matrices 
I< ' in (flTTI and (Tj~5|) have a single non-zero eigenvalue and single non-zero singular value c say 
of degeneracy r, and that c is large. Note that the case r = 1 corresponds to all entries being 
equal to c/N and c/m respectively. Asymptotic formulas obtained in the course of so analyzing 
(jl.ip can then be used to show an analogous effect for (II. 2| . 

Consider first fTTTJ) . We are interested in the p.d.f. for the eigenvalues of H, given the 
eigenvalues {A^} of ij' '. Let this be denoted -P G (A|A^ ^). The Jacobian for the change of 
variables from the entries of the matrix to the eigenvalues and eigenvectors is proportional to 

n \\k-Xjfiuuu) 

l<3<k<N 

(see e.g. |17j ) where (WdU) is the normalized Haar measure for unitary matrices (/3 = 2) and 
real orthogonal matrices (/3 = 1). Hence the sought conditional eigenvalue p.d.f. is proportional 
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to 

(3.1) J] lAfc-A.fe-^^^Af) 2 ) f exp ^ T T(UAU^A^))(uUu) 

l<j<k<N J 

where 

A = diag(A 1 ,...,A JV ), AW =diag(Al 0) ,...,A^). 

Now, with x :— {xj}j—i t ,,, in , y '■= {yj}j=i....,n and (7« 2 \x) denoting the Schur (/3 = 2) 
or zonal (/3 = 1) polynomial indexed by a partition k of n parts, introduce the generalized 
hypergeometric function 

(3.2) rf'tai^Ay" 

It is a fundamental result in the theory of zonal polynomials that (|3.2[) relates to (|3.1[) through 
the evaluation formula (see e.g. [2*5]) 

(3.3) / exp 

and thus the conditional eigenvalue p.d.f. P G (A|A^ ^) is proportional to 

(3.4) J] |A fc -A J f e -^/ 2 )S" 1 ^ + ^ O) ) 2 ) ^ 2//5) (/3A(°);A). 

l<j<fe<JV 

Generalizing (|1.1[) so that it is now proportional to 

(3.5) exp ( - /3Tr(tf - e~ T H^) 2 /2\l - e~ 2T |) 

allows the corresponding eigenvalue probability density function for H , p T say to be written as 
the solution of a Fokker-Planck equation. This is precisely the same Fokker-Planck equation as 
appears in the Dyson Brownian motion model of the Gaussian ensembles (see e.g. [17]). It can 
equivalently be interpreted as the Brownian evolution of a classical gas with potential energy 

1 - 

(3.6) W = - Yl log|A fc - A^-l + 

l<j<k<N j=l 

is given explicitly by 



The Green function solution G^(\\\^) is that solution which satisfies the initial condition 

N 

(3.8) p T (A) _ =I1^- A * (0) ) (A^<---<Al 0) ). 



T = 

1 = 1 



It is given in terms of the generalized hypergeometric function (13. 2|) according to [7] 
(3.9) G^(\\\^) = Ce-e w W 



X6XP V 2{l-P)j^P + {X l > K° V(l-t2)V2'(T- t 2)l/2 



where t :— e T and here and below C is some proportionality constant independent of the 



primary quantities (here {Aj}, {A^}) of the equation. 
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Substituting (|3.9p in (|3 . 3|1 allows the matrix integral to be expressed in terms of the Green 
function. The significance of this for purposes of studying the asymptotics of P G (A|Jv°)) in the 
case that X^ ' is given by 

(3.10) A (0) = ((0) N ~ r , (c) r ) 

is that for r — ► the asymptotic form of G^(A|A'°- ) ) must factorize as 

(3.H) G G ({A i } j=r+1 ,... )iV |(0) JV - I -)GG({A i } j=1> ..., r |( C ) r ) 

since the initial condition (|3.8p then separates the system into two. Furthermore, we know from 
that 

(3.12) G?({A J } i=1 ,.„, n |(cr) -/^re-^^^-^ I** " -M"- 

l<j<fe<n 

Making use of (|3.11[) and (I3.12p in (|3.9p , and making use too of the scaling property 

(3.13) 4 2/f}) (x;ay) = 4 2/P) (ax;y), 

valid for a scalar, allows the following asymptotic expansion to be deduced. 

Proposition 3.1. Let qJq^^ be specified by VS. 2^ . In the limit c — > oo 

n \\k-\jfo4 2/m m((o) N - r Ac) r )) 

l<j<k<N 

(3.14) ^ Ce W^U^-^-f) H \X k -Xjf J] \X k -Xjf. 

l<j<k<r r+l<j<k<N 

For future reference we remark that the next order correction in (|3.14p is to multiply the right 
hand side by the factor (see e.g. [EJ) 

r N 

(3.15) n n 

j=l l=r+l 

which to leading order is a constant. 

Substituting (|3~T4f in (J33J) exhibits the c -> oo form of P G (A|A (0) ). 

Corollary 3.2. The conditional eigenvalue probability density P G (A|<v )), in the case that 
is given by VJ.lCfy with c — > oo, factorizes to be proportional to 

(e-WOEJ-.CA,-)' J] |A fc -A,f)(e-W 2 )^ + ^l J] l^-^f). 

l<j<k<r r+l<j<k<N 

We recognise the first term as the eigenvalue probability density function for the r x r Gaussian 
ensemble centred at A = c, and t/ie second term as the eigenvalue probability density function for 
the (N — r) x (N — r) Gaussian ensemble centred at the origin. 

It is also possible to express the eigenvalue probability density function P L for the general 
variance Wishart matrices (|1.3p in terms of the generalized hypergeometric function qJ 7 ^^ |23j . 
We will revise this point, and then proceed to make use of Proposition ^. ll to deduce a separation 
of eigenvalues in the setting that the covariance matrix E has m — r eigenvalues equal to 1 and r 
eigenvalues equal to 6, for b — * oo. The case r = 1 corresponds to the setting of Proposition ^. 31 

To obtain the generalized hypergeometric function form of the eigenvalue probability density, 
it is most convenient to consider as the input data not the eigenvalues of the covariance matrix 
E, but rather its inverse E _1 . In particular, we are interested in the case that 

(3.16) E- 1 = diag((6) r , (l) m ~ r ) := A^. 
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With the eigenvalues of (|1.3|) denoted {\j}j=i,...,m ancl A : — diag(Ai, . . . , A m ), by making use of 
relevant Jacobians (see e.g. [T7]), P L is seen to be proportional to 

m f R 

(3.17) J] [] |A fc -A,f / exp(-|{/AC/U(°))({/^{/), a = (n-m+l//3)(/3/2)- 1. 

2 = 1 l<j<k<m 

The sought expression involving qJ 7 ^^ now follows by substituting (|3.3|) . thus giving P L as 
being proportional to 

m 

(3.18) l[\r J] |A fc ~A,f ^ 2//3) (A;-(/3/2)A(°)) 

i=l l<j<k<m 

where 

(3.19) A (0) = ((6) r ,(l) m - r ). 

In (|3.18[) . the limit b — > + and thus 6 — > oo in the setting of Proposition ^. 31 does not directly 
correspond to the setting of Proposition 13.11 However, if we first write A = u/b, making use of 
(1 3 . 1 3 [) shows that we seek the asymptotics of 
(3.20) 

4 Vf3 \u- -(/3/2)((ir, (l/b) m - r )) = T,?=i u 3 o jr(V0) iu . _(/J/2)((0)*\ (1/6 - 1)— ')) 

where the equality follows from the matrix integral form of q.Fq 2 ^, (|3.3p . We remark that 
writing A = u/b is well founded because for (|1.3p all eigenvalues are positive, so we expect the 
variables u to be 0(\) in the limit & — > + . Applying Proposition ^. 31 modified to include the next 
order correction term (|3.15p , to (|3.20p and reverting back to the variables A gives the asymptotic 
expansion 

J] |A fc -A,f ^ 2//J) (A;-(/3/2)A(°)) ~ Ce~W l ^U^ 

l<j<k<m 

(3.21) xe -W*)XU +1 x i f[ x WW™-r) Y[ | Afc _ Aj f I] lAfc-A^. 

(=1 l<j<fc<r r+l<i<fe<m 

Substituting this in (|3.19|) we obtain the 6 — > + form of P L . 

Proposition 3.3. XTie eigenvalue probability density P L /or general variance Wishart matrices 
\1.S\) . in the case that is given by A3. 16]) . with b — > + , factorizes to be proportional to 

r rn 
^ X <*+P(m-r)/2 e - {m - bXj -Q | Afe _ A .,/3)( TJ A^"^ [] |A fe - A,f 

J — 1 l<j</c<r j— r+1 r+l<j</c<m 

We recognise the first term as the eigenvalue probability density function for the r x r Laguerre 
ensemble with A i— » &A, and a i— * a+/3(m — r)/2, and thus with support to leading order at 0(l/b), 
and the second term as the eigenvalue probability density function for the (to— r) x (rn— r) Laguerre 
ensemble. 

We turn our attention now to the shifted mean chiral matrices ()1.5|) . in the case that has 
a single non-zero singular value c of degeneracy r. As already noted, the non-zero eigenvalues 
come in ± pairs, and there are n — to zero eigenvalues. Furthermore, the matrix of eigenvectors 
(|2.26p has U nxn ='■ U as the matrix of eigenvectors of YY\ Y := X + 1' ', and V mX m '■= V 
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the matrix of eigenvectors of Y^Y. The Jacobian for the change of variables to the positive 
eigenvalues and eigenvectors is proportional to 

m 

(3.22) II A ' Q+1 II \4 - ^j\ f3 (UUU)(VUv) 

1=1 l<j<k<m 

where a is as in (|3.17|) . 

With X distributed as a Gaussian matrix according to (jl.2p with E = l m , it follows from 
fj3 . 22(1 that the probability density function for the positive eigenvalues of (|1.5p is proportional 
to 



! e /3ReTr([/AVA<°> t ) 



(3.23) JjA? a+1 e-^ A ? [] \4-*jf (tfdU) (VUV), 

1 = 1 l<j<k<m ^ ^ 

where A is the n x m matrix with all entries equal to except those on the diagonal which are 
equal to Ai, . . . A m , and A^ ^ is the mxn matrix with all entries equal to zero except those in the 
first r positions of the diagonal, which are equal to c. We introduce now the further generalized 
hypergeometric function 

(3.24) 



| K |! ci Q) ((l)«) 

where [ffl]« is a certain generalized Pochammer symbol (see e.g. |17j). The matrix integral in 
(13.23|) can be evaluated in terms of o^i according to (see e.g. [17j ) 

(3.25) / (UUU) [ (VUV) e /3RcTr (C /AyA(°)t) = ^/P) { p n/2 . tf. ( /3 /2)2( A (0))2 ) _ 



Hence the probability density function for the positive eigenvalues of (|1.5|) can be written in 
terms of qF^^ , being proportional to 

m 

(3.26) jjA^e-W 2 ^? J] l A ' " A|^o^ (2/ ^(/3n/2; A 2 ; (/3/2) 2 (A(°>) 2 ). 

i=l l<j<k<m 

Suppose now that the distribution of Y := X + X^ is generalized to take on a parameter 
dependent form proportional to 

(3.27) exp ( - /3Tr(Y - e - T X^(Y - e - T X (0 ^/2\l - e - 2r |) 

(cf. (13. 5[> ). It is known (see e.g. [T^]) that the p.d.f. of the positive eigenvalues of Y, with those 
of A"(°) regarded as given, satisfies the Fokker-Planck equation (|3.7j) with 

, » , iv 

(3.28) w = w«* = -x; a 3--5>g a 3- X! logi^i -^i- 

3=1 3=1 l<j<k<N 

Here a' = (2//3)(a + 1/2) and N — m. This relates to the generalized hypergeometric function 
(I3.24p through the fact that the Green function solution of this Fokker-Planck equation, Gf 1 say, 
can be written [7] 

(3.29) Gf(\\\^) = Ce~P w ^ 

where i = e~ 2T as in (|3.9p . and k is independent of A, A' ). 
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We can view (|3.29p as allowing o-7"i in (|3.26p to be rewritten in terms of the Green function. 
But as in ||3~TU|) . with 

(3.30) X® = ((0)™- r , (c) r ) 

the asymptotic form of G^ h (A|A^ ^) must, for r — * 0, factorize as 

(3.31) Gf({Ai}i=r + i 1 „.,m|(0r- r )G? l ({A J } i = 1 ,...,r|(cr). 
Furthermore, we know from [7] that 

n 

(3.32) Gf({A,} j=1 ,...,„|(0n T ~ o CnAf' e-^ 4r ^?- A ? II 1^ - A^, 

3=1 l<j<fc<ri 

while for large c, the fact that with Xj i— > Aj + c (|3.28[l reduces to (|3.6p tells us that 

(3.33) G?({A} i= i,...,r|(cr) ~ G?({A} i=1 ,..., r |(c)-) 

and is thus given by (|3.12[) . This enables us to deduce the following asymptotic expansion, and 
then proceed to deduce the c — > oo form of P ch (A|A'°- ) ). 

Proposition 3.4. Let oJ-"i be specified by {3.2$ . In the limit c — ► oo 

v 

II X ?°' II ^ - A J 2 |^o^ (2//3) (/3n/2; A 2 ; (0/2) W^. (c)-)) 

i=l l<j<k<N 

V 

(3.34) ^ CeW^U^i-^-of) -Q | A 2 fe - A 2 f J] Af Q ' TJ |A 2 .-A 2 ^ 

l<j<fe<r i=r+l r+l<j<fc<V 

Corollary 3.5. XTie conditional eigenvalue probability density P ch (A|A^ ^), in the case that A^ -* 
is given by Iff, fffl)) c — > oo, factorizes to be proportional to 

m 

(e-WmTr^-*? Y[ |A fc -A,f)( f[ \f a 'e-W»$ TJ I**-**!")- 

l<j<fc<r r+1 r+l<j<fc<m 

We recognise the first term as the eigenvalue probability density function for the r x r Gaussian 
ensemble centred at A = c, and f/ie second term as the eigenvalue probability density function for 
the (m — r) x (m — r) chiral ensemble centred at the origin. 

4. Explicit form of the correlations for /3 = 2 

According to (|3.ip and (|3.17p . the matrix integral in (|3.3|) fully determines the eigenvalue 
probability density function for the shifted mean Gaussian and general variance Wishart ensem- 
bles. Furthermore, we know from (|3.23p that the matrix integral in (|3.25p fully determines the 
eigenvalue probability density function for the shifted mean chiral ensemble. In the case (3 = 2 
the matrix integrals in (|3.17p and (I3.25P are over the Haar measure on the unitary group, and 
can be evaluated in terms of determininants (see e.g. [T7]). In fact each of (|3.ip . (|3.17p and 
(|3.23p can then be written as the product of two determinants. From this functional form, using 
the general theory of biorthogonal ensembles [TU] , it is then possible to proceed to compute the 
general n-point correlation functions as an n x n determinant, 

(4.1) P(n)(Ai, ■ ■ A n ) = det[.Kjv(Aj, Xk)]j,k=i,...,n 

for a certain kernel function K^(x,y). In this section the eigenvalue separation phenomenon 
for P = 2, exhibited at the level of the eigenvalue probability density function in the previous 
section, will be analyzed both analytically and numerically in terms of the correlation functions. 
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A 

Figure 1. (Color online) The eigenvalue probability density for the N x N 
GUE for N = 15 with r = 5 shifted eigenvalues, = c = 15, i = 1, . . . , r 
(black solid line). The density has two distinguishable parts, one is similar to 
the eigenvalue probability density for the (N — r) x (N — r) GUE centered at the 
origin (blue dotted line) and the other is similar to the eigenvalue probability 
density for the rxr GUE centered at A = c (red dashed line) . These similarities 
become exact correspondences in the limit c — > oo. 



Consider first the shifted mean Gaussian ensemble, and suppose that XW in d3TXJ) has a single 
non-zero eigenvalue c of degeneracy r, as assumed in Section 3. Define 

(4.2) K%™(x,y) := f^!^! yj gpW^M, 

p; 2V 

where H p (z) denotes the Hermite polynomial. Then we know from [13j that 

r 

(4.3) K N (x,y) = K^(x,y) +J2^ j H^ U) (y), 

3 = 1 

where, with C{ _ 2 c} a simple closed contour encircling zero and -2c, and are so called 
incomplete Hermite functions specified by 



e 



-xz-z 2 /i dz 



/ZOO 1 

According to (|4.ip setting x — y in (|4. 3[) gives the eigenvalue density. This being a function 
of one variable, it is well suited to a graphical representation. Indeed the quantities (|4.2|) and 
(I4.4[) making up (|4.3p are all readily computed numerically, so allowing for particular values of 
TV, r and c the density to be tabulated and then graphed. In regard to (|4.4p . we first make use 
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Figure 2 . (Color online) The onset of eigenvalue separation for the NxN GUE 

for N = 500 with r = 1 shifted eigenvalue, Aj ' = cJ/2 and J = \/2N. Shown is 
the right edge of the eigenvalue probability density for c = (black solid line) , 
c = 1 (blue dotted line), c = 1.2 (red dashed line), and c = 2 (green dot-dashed 
line) . Inset shows the full eigenvalue probability density for the NxN unshifted 
mean GUE with TV = 500. 



of the integral representations of the Hermite polynomial 

e xz-z 2 /4 dz 



H n (x) = 2 n nl 
(4.5) = \/¥e : 



{o} z n+1 2m 

z n e -xz+z 2 /4 _^f_ 

2ni 



to evaluate the contour integrals. Consider for defmiteness f^^(x). First computing the residue 
at z — —2c, then computing the contribution of the singularity at the origin by writing 

i i ^v-n- zP 



2c 2c ^ y ' (2c)p' 

p=Q y ' 



and using the first of the integral formulas in (|4.5p shows 

(4-6) rW(*) = (-i)^(^ r -X; 



2cx ~ c2 — 2 1 H N ^ p (x) 



(2c) N ~ 1 ^ (2c)p+! 2 n - 2 -p(N -2-p)\ 

An example of a numerical calculation illustrating Corollary 13. 21 is given in Figure [T] Another 
application is to use a numerical evaluation of (|4.3[) to illustrate the onset of the eigenvalue 
separation for large N as specified in Proposition 12.21 This we do in Figure [H For these 
numerical calculations, we evaluated the integrals in (4.4j) for arbitrary j obtaining expressions 
involving sums of Hermite polynomials as in (|4.6|) and used them in (|4.3|) . 

For general x, y we can show analytically that the large c form of pi n \ is consistent with 
Corollary 13. 21 by establishing the following asymptotic result. 



16 



K.E. BASSLER, P.J. FORRESTER, AND N.E. FRANKEL 



Proposition 4.1. For c — > oo 

r 

(4.7) J2r^(x)T^(y) ~ e 2c ^-^^ UE (^-c,j/-c). 

3=1 

Proof. We see from the first formula in (14. 4|) that 

e -:rz-2 2 /4 ^ z g2ca;-c 2 /■ — (x-c)z-z 2 /4 ^ z 



f(j')( 



(_l)J-l e 2cx-c 2 Hj-iix-c) 



zi 2m 

{-2c} 



(-2c)^- r 2i~ 1 {j-l)\ 

where the second equality follows from the first integral formula in (|4.5p . For the second integral 
formula in (|4.4|) , use of the second integral form of the Hermite polynomials in (|4.5| shows 

rW)(y) ~ (-l)^ 1 (-2c) jv -'-e^ 2 g3 ' l( ) l~ C) . 

c— >oo ^/7f 

Substituting in the LHS of (|4"77|) and comparing with (gjg) gives the RHS of P~7) . □ 

According to ()4.3p and l|4.7jl . for c — > oo the correlation kernel Kn is the sum of two terms, one 
with support in the neighbourhood of the origin, and the other with support in the neighbourhood 
of x, y — c. Furthermore, these two terms are the correlation kernel for the (N — r) x [N — r) 
GUE and the r x r GUE, the latter shifted to be centred about A = c (note that the factor 
e 2e(x-y) m (|4 joes note effect the determinant (14. 

In the case of general variance complex Wishart matrices (II .3[) with £ _1 specified by (|3.16p , 
the correlations are given by (|4.1|) with N = m and 

r 

(4.8) K m (x,y)=KZ? r (x,y)+ 1 £ i A®(x)A®(y) 

i=i 

(see |13]V Here a = n — m, and with L^{x) denoting the Laguerre polynomial 

(4.9) K«(x,y) = y a e-yJ2 ij ^^L;(x)L;(y) 

while AW and A^ have been termed incomplete multiple Laguerre functions and are specified 
by 

AW J e-"(l + ,r dg 

•/c { o,(S- 1 )} ZTO " r ( 2 -^- 1 )) i27ri 

(4 10) AW (x) - I ( W -(6-l)r^ du. 



Making use of the integral representation for the Laguerre polynomials 

w n+1 v ' 2m 



(4.H) £J(i)= ^r(i + ^r +Q ^- 



-{0} 

allows AW (x) to be expressed in terms of these polynomials. For example, by considering 
separately the neighbourhoods of zero and b — 1 as in the derivation of (|4.6p , we can use (|4.1ip 
to show 

(6-1) — r ( 5_ 1)P+ 1 • 
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FIGURE 3. (Color online) The eigenvalue probability density for the m x m 
variance LUE or variance complex Wishart matrix ensemble for m = 10 and 
a = 1/2 with r — 3 shifted eigenvalues, A- ^ = b = 0.05, i = 1, . . . ,r (black 
solid line). The density has two distinguishable parts, one part, shown in (a), is 
similar to the eigenvalue probability density for the unshifted (m — r) x (m — r) 
variance LUE with unchanged a (blue dotted line) and the other part, shown in 
(6), is similar to the eigenvalue probability density for the r x r variance LUE 
with A i— > bX and a i— > a + (m — r) (thick red dashed line). These similarities 
become exact correspondences in the limit b — > + . Note that in (a) the two 
lines (black solid and blue dotted) are indistinguishable and that the shifted part 
of the density is off-scale to the right and not visible in the figure, while in (b) 
a difference between the two lines (solid black and thick red dashed) is slightly 
apparent. The large sharp peak at small values of A in (b) is the unshifted part 
the density shown in (a). 
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Figure 4. (Color online) The onset of eigenvalue separation for the m x m 
variance LUE or variance complex Wishart matrix ensemble for m — 500 and 
a = 1/2 with r = 1 shifted eigenvalue, A± = bJ and J = 4m. Shown is the right 
edge of the eigenvalue probability density for b = (black solid line), b = 0.5 
(blue dotted line), b = 0.45 (red dashed line), and b — 0.275 (green dot-dashed 
line). Inset shows the full eigenvalue probability density for the m x m unshiftcd 
mean variance LUE with m — 500. Note the logarithmic scale of the vertical 
axis of the inset graph. 



To relate to La guerre polynomials requires the identity 

n\(-x)- a L- a (x) = (n - a)lL a n _ a (-x), a e Z, 
which when used in (|4.11[) implies the integral representation 

e*x- a (n + a)l f e™ n dw 



c { _ 1} fa + l) n+a+1 



From this we obtain, for example, 



(m — 1)! 

With the integrals (|4. 10|) thus made explicit for arbitrary j, (j4. 8[) in the case x = y, and for 
particular m, a and r can readily be tabulated and graphed. An example illustrating Proposition 
3.31 is given in Figure [3] However, unlike the situation with (14. 3p . the structure (|4.8|) is not well 
suited to exhibit the eigenvalue separation effect of Proposition 13.31 in a single plot. This is due 
to the dependence on a + r rather than a in the first term on the RHS. We can also make use of a 
numerical evaluation of (|4.8p in the case x = y to illustrate the onset of the eigenvalue separation 
for large m and b near 1/2, as predicted by Proposition ^. 31 This is done in Figure [U 

It remains to consider the shifted mean complex chiral matrices (jl.5|) in the case that X^ ^ 
has a single non-zero singular value c of degeneracy r as in Section 3. According to [14] the 



EIGENVALUE SEPARATION IN SOME RANDOM MATRIX MODELS 19 




Figure 5. (Color online) The eigenvalue probability density for the m x m chiral 

matrix ensemble for m — 15 and a = 4 with r = 5 shifted eigenvalues, A^ = 
c = 15, i = 1, . . . , r (black solid line). The density has two distinguishable parts, 
one part is similar to the eigenvalue probability density for the (m — r) x (m — r) 
unshifted mean chiral matrix ensemble with the same a (blue dotted line) and 
the other part is similar to the eigenvalue probability density for the r x r 
GUE centered at A = c (red dashed line). These similarities become exact 
correspondences in the limit c — * oo. 



correlations are specified by 

n 

(4.12) P(i)(xi, ■ ■ ■ , xi) = 2 l Y[ x v det[if m (a;|,4)kfc=i,...,i 
where 

r 

(4.13) K m (x,y) = K^_ r (x,y) + Y / Mx)q l (y). 

i=l 

In !j4l5J) K^_ r is specified by 0~U), while 

X /•oo 

(4.14) p k (x) = rf — u' n+a - r {u + c 2 ) k - 1 e- u F 1 (a + l--xu)du 

r(a + l) J 



(4.15) q k (x) 



e v oFi(a + I; —xv) dv 
T(a + 1) J e _j _ o3 w m - r (-y + c 2 ) fc 2?ri' 



The functions Pk(x), qk(x) can be expressed in terms of Laguerre polynomials. To see this, 
for Pk(x) we require the integral formula 



(4.16) L°(s)= — / u^e^o^^+l;-^)^, 

n!r(a + 1) 7 

while for gfe(x) we require 

L*(a;) _ 1 f e w F 1 (k + l;-xw) dw 



T(n + k + l) r(fc + l)7 C{o} w n+1 2m 
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Figure 6. (Color online) The onset of eigenvalue separation for the to x to 
chiral matrix ensemble for to = 500 and a = 2 with r = 1 shifted eigenvalue, 
= cJ/2 and J = 2-^/m. Shown is the right edge of the eigenvalue probability 
density for c — (black solid line), c = 1 (blue dotted line), c = 1.2 (red 
dashed line), and c = 2 (green dot-dashed line). Inset shows the full eigenvalue 
probability density for the to x m unshifted mean chiral matrix ensemble with 
m = 500 and a = 2. 



Once Pk{x), qk{x) have been expressed in terms of Laguerre polynomials, the kernel (|4.13p is 
for x — y and particular to, a, r readily computed numerically. A plot illustrating the eigenvalue 
separation property Corollary 13.51 is given in Figure [5j A further plot illustrating the onset of 
the large to separation of a single eigenvalue, for appropriate c as specified in Proposition ^. 51 is 
given in Figure El 

The correlation (14.12p for general / can be shown analytically to be consistent with Corollary 
3.51 For this the c — > oo behaviour of the summation (14.13|) is required. 



Proposition 4.2. For large c 



(4.18) 



qk(x) 



(-1) 



N-r+k-l 



-H k -i(Vx- c) 



c 2(N-r)+k 0f2 fc (fc-l)! 

while for large c with y/x — c fixed 

(4.19) pk(x) ~ {-if-^c^-^^Hk-^- c) 

Hence for large c 



(4.20) 



2c 



-K 



GUE 



{y/x- c, A /y- c). 



Proof. We require the formula 
(4.21) I a (2^) 



r a/2 



r(a + l) 



Fi(a + l;z) 
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and the asymptotic expansion 

(4.22) I a (x) ~ e \ 

x^Q (2lTX) L ' Z 

Consider first {x) ■ For c large we see that the main contribution comes from the singularity at 
v = — c? . Writing v i— > — c 2 + v, making use of (|4.21|) and (|4.22p . then expanding the exponent 
to second order in v while keeping only the leading term in the rest of the integrand shows 

c -a x u/2 e -(^x-c) 2 y f e -(\/x-c)v-v 2 /4 dv 



2y/^(Vx^) 1 / 2 (-C 2 ) N - r C k - 1 J C{0} V k 2%l 

The contour integral can be evaluated in terms of a Hermite polynomial according to the first 
formula in (14. 5j) . and furthermore x can be replaced by c 2 in the prefactors not involving the 
difference y/x — c. This gives (|4.18[) . 

Consider now pk(x). Making use of (|4.21|) and (|4.22|) . then expanding the exponent to second 
order about the stationary point at u = — x and deforming the contour into the direction of 
steepest descents shows 

2(N-r)+k-l foe 

p k (x) ~ (-l) N ~ r — j= / {iu ~ 2(v^ - c)f- x e- u / 4 du. 



2v^ 

Recalling now the second integral formula for the Hermite polynomials in (|4.5[) gives (14. 19|) . 
Substituting P~T5|) and (f4~T^ in the LHS of ([430]) . and recalling (O gives the RHS of fl~20l) . 

□ 

Substituting (f4~^0|) in (|4"T2"|) we see that for c — > oo the correlation kernel separates into two 
parts, and that the resulting correlations are those corresponding to the two ensembles identified 
in Corollary [ 



5. A FURTHER ASYMPTOTIC LIMIT FOR THE GENERAL VARIANCE LAGUERRE ENSEMBLE 

We see from ()3.18j) that the general variance Wishart matrices (|1.3j) . in the case that is 
specified by (|3.16p . is well suited to asymptotic analysis of the b — > oo limit. In the setting of 
Proposition 12.31 this corresponds to b — > 0, so the analysis leading to the eigenvalue separation 
therein has no bearing to this limit. Starting with ()4.8|) . and with m — * oo, the scaled correlation 
functions corresponding to this limit were analyzed in (13] . The scaling was specified by setting 
Xi Xi/m and b — mc, c fixed. Here we will consider the limit b — > oo with m fixed. 

Proposition 5.1. The eigenvalue probability density P L for general variance Wishart matrices 
in the case that S _1 is given by 113.16]) . with b — > oo, factorizes to be proportional to 



r m 

n>?^ n i^-A,r)( n n i a *-v 



j = l l<j<k<r j=r+l r+l<j<k<m 

We recognise the first term as the eigenvalue probability density function for the r x r Laguerre 
ensemble with A i— > bX and thus with support to leading order at 0(1/6), and the second term 
as the eigenvalue probability density function for the (m — r) x (m — r) Laguerre ensemble with 
ot i ► ol ~\~ r . 

Proof. This is an immediate consequence of the asymptotic formula (|3.20[) , modified to include 
the correction term (|3.15|) . □ 



In the case (3 = 2 the eigenvalue separation property can be exhibited numerically using (|4.8|) . 
For this we proceed as in the lead up to Figure [31 making use of evaluations in terms of Laguerre 
polynomials of the quantities J4.10|) . An example is given in Figure [7] 
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Figure 7. (Color online) The eigenvalue probability density for the m x m 
variance LUE or variance complex Wishart matrix ensemble for m = 20 and 
a = 3 with r = 5 shifted eigenvalues, X^' =b= 100, i = 1, . . . , r (black solid 
line). The density has two distinguishable parts, one part, shown in (a), is 
similar to the eigenvalue probability density for the unshifted (m — r) x (m — r) 
variance LUE with ana+r (blue dotted line) and the other part, shown in (6), 
is similar to the eigenvalue probability density for the r x r variance LUE with 
) h M and unchanged a (thick red dashed line). These similarities become 
exact correspondences in the limit b — > oo. Note that in (a) the two lines (black 
solid and blue dotted) are almost indistinguishable and that the large sharp peak 
at small values of A, which does not correspond to the blue dotted line, is the 
shifted part the density shown in (b). In (b) the unshifted part of the density is 
off-scale to the right and not visible in the figure. 
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6. Concluding remarks 

The prime motivation for our work was to advance the understanding of the precise effect 
shifting the mean in the Gaussian probability distribution for the elements ([l.ip has on the 
eigenvalue distribution. We sought to compare this to analogous effects for general variance 
Wishart matrices, and shifted mean chiral ensembles. We also had in mind the application of 
our work in the fields of biological webs, food chains, plant and animal ecology, and networks, 
neural and otherwise. 

During the period of the 1950s through the 1970s there was an intense activity by ecologists in 
studying the role played by biodiversity on ecological systems. In particular, how the number of 
species of a community affected its stability. There was quite conflicting evidence as to whether 
complexity increased or decreased or had no effect on the stability of an ecosystem. 

The coup de grace came with the classic work of May [301 1311 132] from purely a theoretical 
approach. This work has inspired literally a vast amount of literature on the complexity-stability 
issue in these and other arenas. May used the notions and technical details from the mathematics 
of random matrix theory in formulating his simple as he called them mathematical models with 
their predication that increasing complexity decreases the stability of the system. 

This pioneering work has motivated in the past 35 years and will continue to do so into the 
future the debate on the complexity- stability issue. In particular, it has been evidenced that 
under certain specific conditions in the system that complexity can either increase or decrease 
the system's stability. In an interesting toy model [21] the authors have studied an ecological 
system where they use a small matrix where the elements are taken from a random distribution 
of interaction coefficients. They found that the non zero mean of the distribution can lead to 
either an increase or decrease in the stability of the system. (See also some very early work of 
1971 [33J discussed by May [3JE2].) 

Although May's work was strongly influenced by random matrix theory, there has been little 
application of the mathematics of random matrix theory to these vitally interesting areas and 
issues. Specifically, it is the distribution of the largest eigenvalue(s) of a random matrix that 
determines its stability. As we have shown in this paper, a non zero mean in the probability 
distribution for the elements of a random matrix severely affects its eigenvalue distribution and 
hence its stability. 

Very recent works [HI [37l [21 E] studying the eigenvalues of the real Ginibre ensemble [21] 
have strong potential bearing on the further application to biological webs and neural networks 
[361 [TB, 39J. The generalization to the case of complex eigenvalues is very important and we plan 
to study this in our future work. 
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